Minimizing the linewidth of the Flux-Flow Oscillator 
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For the first time the linewidth of Flux-Flow Oscillator has been calculated by direct computer 
simulation of the sine-Gordon equation with noise. Nearly perfect agreement of the numerical 
results with the formula derived in [Phys. Rev. B, 65, 054504 (2002)] has been achieved. It has 
been demonstrated that for homogeneous bias current distribution the linewidth actually does not 
depend on the junction length for practically interesting parameters range. Depending on the length 
of the unbiased tail, the power may be maximized and the linewidth may be minimized in a broad 
range of bias currents. The linewidth can be decreased further by 1.5 times by proper load matching. 
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During the last decade the flux- flow oscillator (FFO), 
based on a viscous flow of magnetic flux quanta in a long 
J osephson tunnel junction ( JT J) [l[ , has been considered 
as the most promising local oscillator in superconducting 
spectrometers 0] for space-born radio astronomy and at- 
mosphere monitoring due to its wide operational band- 
width and easy broadband tunability. However, the spec- 
tral linewidth of the emitted radiation of the free-running 
FFO is rather large, that complicates the phase locking. 
Typically, the free-running linewidth is 2-20 MHz for an 
Nb-AlOx-Nb FFO in the 400-700GHz frequency range. 
For spectral applications it is of crucial importance to re- 
duce the FFO linewidth, to make it more homogeneous in 
all working frequency range, and to increase the emitted 
power to improve the signal-to-noise ratio. 

The dynamical and fluctuational properties of the 
FFO have been investigated in [l|]-[2(3], m particular, the 
linewidth has been studied both experimentally 
and theoretically 14|- [l8j | . H owever, the only formula for 



the linewidth, derived in 17|, that takes into account not 



only differential resistance over bias current, but also dif- 
ferential resistance over magnetic field, has been proven 
as the most adequately describing experimental results, 
see 11 1 . But, the nature of conversion of bias cur- 
rent fluctuations to the magnetic held fluctuations is still 
unclear, and the conversion factor is not known exactly. 
Also, the dependence of the linewidth on the bias cur- 
rent profile and certain parameters, such as RC-load, has 
not been systematically studied yet neither theoretically 
nor experimentally. The aim of the present paper is to 
study the FFO linewidth by direct computer simulation 
of the sine-Gordon equation with noise, to compare the 
obtained results with the formula of [17| and to make 
certain optimizations of FFO design (varying bias cur- 
rent profile and the RC-load) in order to minimize the 
linewidth and to increase the emitted power. 

For several decades the sine-Gordon model has been 
the most adequate model for the long Josephson junc- 
tion, giving a good qualitative description of its basic 
properties: 

4> tt + a(j> t - (j> xx = f34> xxt + T)(x) - sin(0) + r)f(x, t), (1) 



where indices t and x denote temporal and spatial deriva- 
tives, respectively. Space and time are normalized to 
the Josephson penetration length Aj and to the in- 
verse plasma frequency w" 1 , respectively, a = uj p /uj c is 
the damping parameter, where u> p = y / '2eI c /hC , u> c = 
2eI c Rjq /%, I c is the critical current, C is the JTJ capac- 
itance, and i?jv is the normal state resistance, (3 is the 
surface loss parameter, rj{x) is the dc overlap bias cur- 
rent density, normalized to the critical current density J c , 
and rjf(x,t) is the fluctuational current density. In the 
case where the fluctuations are treated as white Gaussian 
noise with zero mean, and the critical current density 
is fixed, its correlation function is: (if(x,t)if(x' ,t')) = 
2a r yS(x — x')S(t — t'), where 7 = It/(J c ^j) is the dimen- 
sionless noise intensity [2(J, It = 2ekT/h is the thermal 
current, e is the electron charge, h is the Planck constant, 
k is the Boltzmann constant and T is the temperature. 

The boundary conditions, that simulate simple RC- 
loads, see Ref.s [3[ and [19j, have the form: 

0(0, t) x + r L c L 0(O, t) xt - c L 0(O, t) tt + (2) 
0r R c R (j>{Q, t) xtt + /30(O, t) xt = r - Ar, 
<j>(L, t) x + r R c R 4>{L, t) xt + c R (j)(L, t) tt + (3) 
Pr R c R( j){L, t) xtt + p<f>(L, t) xt = r + Ar. 

Here T is the normalized magnetic field, Ar = 0.05r 
is a small magnetic field difference, see [19j, and L is 
the dimensionless length of JTJ. The dimensionless ca- 
pacitances and resistances, cl,r and rj,^, are the FFO 
RC-load placed at the left (output) and at the right (in- 
put) ends, respectively. It should be noted that, following 
Ref. [l8j], if both overlap i] ov = (l/L) rj(x)dx and in- 
line r]i n = 2Ar/L components of the current are present, 
the total current, rj t , with respect to which all current- 
voltage characteristics will be computed, is the sum of 
overlap and inline components: rjt — t]ov + Tjin- 

In Ref. [3] on the basis of Eq. JT]) without noise and 
boundary conditions ^ , © , the investigation of current- 
voltage characteristics of FFO has been performed. For 
the bias current profile, depicted in the inset of Fig. [l]by 
curve with crosses, good qualitative agreement with ex- 
perimental IVCs has been achieved. Due to experimental 
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motivation it was assumed, that the current profile was 
parabolic (with the curvature a = 0.005) between the 
left and the right boundaries of bias electrode Xq and X\ 
(0 < xq < x\ < L), and drops down exponentially in the 
unbiased tails x < xq, x > x\ with the decay factor p: 
cxp(-px) (with p = 0.13 in Fig. 1). The decay factor 
and the parabolic curvature were used as fitting param- 
eters when the comparison with the experimental IVCs 
was done. 




FIG. 1: The distribution of overlap component of bias cur- 
rent r)(x). Short-dashed line - homogeneous distribution; solid 
curve - xi = 30; curve with rectangles - xi = 25; dot-dashed 
curve - x\ = 20; long-dashed curve - X\ = 10. Inset: curve 
with circles - xo — 11, x\ = 25.5, a — 0; curve with crosses - 
x = 11, xi = 25.5, a = 0.005. 



The key question arises: if one will change the lengths 
of unbiased tails, how it will affect the emitted power 
and the linewidth? In Ref. [l[ it was suggested to use 
the unbiased tail to decrease the differential resistance r<z 
that might reduce the linewidth, if the formula for the 
linewidth of short JTJ would work for FFO (here and 
below the linewidth is defined as full width, half power): 



A/ s = 2 ai r 2 JL. 



(4) 



Later, it was found experimentally [9(, that even for small 
Td the FFO linewidth is almost one order of magnitude 
larger than predicted by (HJ) . In Ref. (l7[ , the formula for 
the FFO linewidth, that takes into account not only con- 
ventional differential resistance over bias current, but also 
differential resistance over mag netic field r% L = Ldv/dT 
(control line current) was derived: 
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and demonstrated good agreement with experiment 
[Tol | , 11 1 . In Ref. it has been demonstrated that by 
the choice of the bias current profile, the radiation can 
be either enhanced or suppressed, and it is desirable to 
supply more bias current at the radiating end, than at 



the input end. Therefore, in the frame of the present pa- 
per we shift the current profile to the left, xq = 0, and 
vary the length of unbiased tail, which is located at the 
right end of JTJ. Also, to avoid problems with scaling of 
parabolic curvature (it is not clear what to keep constant, 
the curvature or the depth of the "well"), let us set the 
current profile to be constant between xq and X\. 




FIG. 2: Radiated power versus total current, computed for 
rj(x), presented in Fig. [l] short-dashed line - homogeneous 
distribution; solid curve - xi = 30; curve with rectangles - 
xi = 25; dot-dashed curve - xi = 20; long-dashed curve - 
xi — 10; curve with circles - a;o = 11, x\ — 25.5, a — 0; 
curve with crosses - xo — 11, xi = 25.5, a = 0.005. Inset: dc 
current- voltage characteristics computed for r/(x), presented 
in Fig. [T] the notations are the same as for power. 

The power at RC-load at the radiating end x = for 
different bias current profiles depicted in Fig. [H is pre- 
sented in FigUl The power is computed in accordance 
with 0]. The implicit difference scheme, used to solve 
Eq. (|T|) with noise, has been successfully tested in [2(| 
when the mean escape times from zero voltage state were 
investigated. The parameters are the following: L = 40, 
a = 0.033, (3 = 0.035, c L = c R = 100, r L = 2, r R = 100, 
r = 3.6, and 7 = 0.1. From Fig. [5] one can see, that 
for the case of two unbiased tails xq — 11, x\ — 25.5 the 
power is minimal and almost order of magnitude smaller 
than for the current profile with the one unbiased tail 
xi = 25, which gives maximal power among all consid- 
ered current profiles. In the inset of Fig. [5] the current- 
voltage characteristics for the same current profiles are 
given for comparison: it is seen, that the flux-flow steps 
have largest height also for x\ — 25. The height of IVCs 
for both profiles with xq = 11, x\ = 25.5 and a — 0.005, 
a = have close values to each other and are compa- 
rable to X\ — 25. So, it is desirable to give larger bias 
current at the radiating end to get larger emitted power. 
It should be noted that for both profiles with two un- 
biased tails, the pow er versus bias current has minima 
as in experiment [13j. However, the investigation of this 
phenomenon is out of scope of the present paper and will 
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be presented elsewhere. 
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FIG. 3: The calculated spectral density for x\ = 25, r/o = 0.3, 
7 = 0.1 - diamonds; Lorentzian approximation - solid curve. 
Inset: the same, enlarged around the spike. 

The power spectral density of FFO was computed 
as Fourier transform of the correlation function of the 
second kind $[t] = 7^— f Q av (v Q (t)v (t + r)) dt, where 
v (t) — dip(t, 0)/dt is the voltage at the RC-load (x = 0) 
and T av is the averaging time. There are two general 
restrictions, that complicates the calculation of the spec- 
tral density and the linewidth: on one hand, the time 
step should be small enough to resolve oscillations, and 
the averaging time T av should be rather large to resolve 
fine spectral spikes. Due to these restrictions the noise 
intensity was chosen 7 = 0.1. Nevertheless, this is the 
same limit of low noise intensity as in experiments, since 
IVCs are almost unaffected by noise, the spectral spikes 
are narrow, and the linewidth perfectly scales proportion- 
ally to the noise intensity, see below. For the linewidth 
calculations the following parameters were used: spatial 
step Ax = 0.05, temporal step At w 0.1, T av = 8000. 

In Fig. [3] the power spectral density of FFO is pre- 
sented. As one can see, the emitted signal is nearly si- 
nusoidal, in agreement with Q and experimental results: 
the power contained in the second and third harmonics 
is much lower than in the main one. Also, the spectral 
peak is perfectly Lorentzian in more than two orders of 
magnitude interval. 

In Fig. 2] the FFO linewidth versus differential resis- 
tance r<i for junctions of L = 40 and 7 = 0.1 is presented 
for the case of homogeneous bias current distribution and 
x\ = 30. It is seen, that in both cases formula (O is in 
good agreement with numerical results, while formula (j4|) 
and formula from (in formula Q a must be substi- 
tuted by inverse static resistance r) t /v) significantly un- 
derestimate the linewidth. 

The appearance of the conversion of bias current fluc- 
tuations to magnetic field fluctuations may be explained 
in the following way: boundary conditions |J5J|, ((3J) of 
Eq. fl} depend on the phase, which fluctuates, since it 



4H 




FIG. 4: FFO linewidth versus differential resistance ra for L = 
40 and 7 = 0.1. Empty triangles and crosses - simulations and 
theory ([5| for homogeneous bias current distribution, empty 
circles and crests - simulations and theory for the case with 
unbiased tail Xi = 30, filled diamonds and rectangles - theory 
of [Hj], solid line - theory @. 
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FIG. 5: FFO linewidth versus total current, computed for 
r](x), presented in Fig. [TJ 7 = 0.1: short-dashed curve - 
homogeneous distribution; solid curve - xi = 30; rectangles 

- Xi = 25; circles - xi = 20; long-dashed curve - Xi = 10; 
curve with crosses - xo = 11, xi = 25.5, a — 0.005. Inset: 
FFO linewidth for junctions with different lengths: L = 20 

- diamonds, L — 40 - crosses, L = 60 - circles, L = 80 - 
triangles, solid curve - theory for L = 60, all these curves for 
7 = 0.1; rectangles - simulations for L — 40, 7 = 0.05. 

is governed by Eq. |T]) with noise. It is quite difficult 
to calculate the noise conversion factor a, since it should 
be done self-consistently. In the following, to see how 
the linewidth behaves in all working range, let us plot it 
versus total bias current r\t- 

Let us analyze the dependence of the linewidth on the 
length of unbiased tail. From Fig. [5] it is seen, that min- 
imal value of A/ is reached for several cases, including 
the homogeneous one. However, for the unbiased tail 
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x\ = 30, the linewidth is nearly constant (and minimal) 
in a maximal range of bias currents. Therefore, the un- 
biased tail of 1/4 of junction's length gives the nearly 
maximal power and nearly minimal linewidth in broad- 
est range of bias current, and can be recommended for 
spectral applications. 

Let us investigate the dependence of the linewidth ver- 
sus junction's length L for homogeneous bias current dis- 
tribution. The corresponding curves are presented in the 
inset of Fig. [5l It is seen, that increase of L does not 
help to decrease the linewidth, the curves nearly coincide. 
The lowest curve is computed for noise intensity 7 = 0.05 
and the linewidth is two times smaller than for 7 = 0.1. 
Therefore, the noise intensity 7 = 0.1 is indeed rather 
low, and in the following one can get good estimate for 
experimental parameters by multiplication of the com- 
puted curve on the numerical factor 7 e /7, where 7 e is 
the noise intensity, corresponding to the experiment. 
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FIG. 6: FFO linewidth versus total current computed for 
different values of load resistances and noise intensity, x\ — 
30: solid curve - tl = 2, m — 100; triangles - tl = tr = 8; 
crosses - Tl = 2, tr — 8; short-dashed curve - rt = tr = 4; 
curve with diamonds - tl — tr = 2; circles - ri, = t\r = 1; all 
these curves for 7 = 0.1; long-dashed curve - tl = 2, tr — 100 
for 7 = 0.05. Inset: the coefficient a versus tr, rectangles - 
tl — tr; circles - tl — 2; solid curve - fitting by \fr~R. 

Normally, FFO radiating end x = is well matched 
to the external environment, while the opposite end is 
strongly mismatched, as it was modeled in [19( and in 
the present paper. It is interesting to analyze, how the 
linewidth will change, if the FFO is better matched at 
both ends. The results of this analysis are presented in 
Fig. E] for x\ = 30. It is seen, that improved matching 
at the opposite end can decrease the linewidth by 1.5 
times. It is important to note, that equal matching at 
both ends gives smaller linewidth than perfect matching 
at the radiating end and bad at the opposite one, e.g., 
compare the curves for rt = 2, rji = 100 and rt = 
rji — 4 and note, that the curves for rt = Tr = 8 and 



Tl = 2, m = 8 nearly coincide. Finally, it is important 
to mention, that the noise conversion factor a perfectly 
scales as *Jt~r, both for rt = 2, and for rt = tr. 

The linewidth of Flux-Flow Oscillator has been cal- 
culated by direct computer simulation of the modified 
sine-Gordon equation with noise which takes into account 
surface losses and RC load. Nearly perfect agreement of 
the computer simulation results with the formula ([5|) has 
been achieved. It has been demonstrated that for homo- 
geneous bias current distribution the linewidth actually 
does not depend on the junction length for practically 
interesting parameters range. Depending on the length 
of the unbiased tail, the power may be maximized and 
the linewidth may be minimized in a broad range of bias 
currents. The linewidth can be decreased further by 1.5 
times by proper load matching. 
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